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Point to point vacuum correlators containing diquarks in the color anti-triplet representation are 
computed both in the quenched approximation and dynamical overlap simulations with two fla- 
vors. The scalar, pseudoscalar and axial vector diquarks are combined with light quarks to form 
color singlets. The scalar ("good") diquark channel shows a stronger attraction than the axial 
vector ("bad") channel in the quenched data set. The pseudoscalar diquark channel shows a finite 
volume zero mode artifact: the correlator becomes negative at large distance when the quark mass 
is small. By separating configurations without zero modes from those with zero modes, we found 
that the zero modes have an important contribution in both the attraction in the scalar channel and 
the repulsion in the pseudoscalar channel. In the axial vector diquark channel, we did not find 
apparent zero mode effects. 
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1. Introduction 



It is possible that bound states of two quarks, diquarks, may be a component of the structure 
of hadrons. For a review of their phenomenology, see Although the constituent quark model 
has enjoyed a rather successful description of a huge body of hadronic states, it cannot explain all 
aspects of hadron phenomenology. One example is the paucity of exotics which are states beyond 
1717^ baryons and qq mesons. As a special case: why do not a proton and a neutron in close contact 
merge into a single state? Recently, the possible existence of exotic particles ^ |3|] such as the 
pentaquark [Q] has motivated research, e.g. ^, ^ |8[ ^, to understand the role of diquarks in 
QCD. 

In this work, we study diquarks by investigating the properties of point to point vacuum baryon 
correlators in quenched and two flavor dynamical overlap simulations. We calculate several current 
correlators in which two quarks are combined into scalar, pseudoscalar and axial vector diquarks 
in the color anti-triplet representation. Each correlator is normalized by its free field theory version 
so that we can see if there is attraction or repulsion in a specific channel. For each channel in 
the quenched data set, we separate the configurations with zero modes from those without zero 
modes. By comparing the correlators from the two groups, we can try to detect zero mode effects. 
The dynamical simulation is done with fixed topological charge 2 = or ±1. i.e., there is no 
zero mode or one zero mode with positive (negative) chirality in each data set. By comparing the 
correlators from different Q sectors, we can again examine the zero mode effects. 



2. Methodology 

Diquarks are not color singlets. To investigate their properties on the lattice, a diquark can be 
combined with a static (infinite heavy) quark to form a colorless state as was done in Refs. [^, |8]]. 



In Ref. []11[], gauge dependent diquark correlation functions were calculated in a fixed gauge, the 
Landau gauge, on the lattice. Whatever a diquark is, the environment it feels in a baryon with one 
heavy quark is different from that in a light baryon with three light quarks. Here we combine a light 
quark with a diquark with a specific quantum number in the color anti-triplet representation to get a 
color singlet and compute the point to point correlation function. By doing this, we cannot extract 
the mass splitting between spin and spin 1 diquark states or the size of a diquark state as were 
done in Refs. [0, |, ^ since the interaction between a diquark and a light quark also depends on 
the spin of the diquark. However, by comparing baryon correlators containing different diquarks, 
we can see which diquark is favored and which is not. We can also investigate the effects of the 
zero modes. Excess zero modes of Dirac operators are artifacts in quenched simulations. 

The currents and correlation functions we considered are collected in Table [I[ Here C is the 
charge-conjugation operator. We choose not to consider diquarks in the color sextet representation 
because they have much larger color electrostatic energy and thus are not favored phenomenolog- 



ically []12[]. In perturbative QCD, one-gluon exchange leads to attraction [ ]13| , |14| ] between two 
quarks in the color anti-triplet representation with 7^ = 0+. Instanton interactions [15, 16, 17] is 
also attractive in the scalar diquark channel. Because it is thought to be attractive, this channel is 
called the "good" diquark. We use the current in Table ^ to investigate this channel. In quenched 



simulations [18], the point to point scalar meson correlator was found to be negative, a quenching 
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Table 1: Currents and correlation functions 



(diquark) 


Color 


Current 


Correlator R(x) 


0+ 


3 




iTr[(a|r75(x)75(o)|a)xv7v] 





3 


j' = eahcWCd^]u' 


^Tr[(a|r/(x)/^(0)|n)xv7v] 


1+ 


3 




iTr[(n|r73(x)/3(0)|a)xv7v] 



artifact. We are curious to see if we will see similar artifact in the pseudoscalar diquark channel by 
using 7^. (This effect is predicted by instanton liquid models.) The 1+ diquark is called the "bad" 
diquark in the literature, since all models suggest that it is heavier than the scalar diquark. 

The two point correlator for a current / is defined as (n|r/(x)/(0)[n), where |n) is the vac- 
uum and r is the time order operator. For free massless quarks, the quark propagator in coordinate 
space takes the form (in Euclidean space) 

(0|r^(x)^(0)|0) = 2^^g^. (2.1) 

Here \i is summed over. Then for the current in Table ^ we have 

(0|r75(x)75(0)|0) = -^^. (2.2) 

Similarly, we can get the free correlators for other currents in Table They are the same as the 
result in Eq.(^^) except for a sign flip for (for 7^, X3 = is needed to get a same result). As was 



done in Ref. [ 19, 2C ], it is convenient to multiply the correlators with Xv7v and take the trace in the 



Dirac indices to get a number for each of the currents. For example, from Eq.(I^), we find Rq{x) = 
^Tr[(0|r7^(x)/^(0)|0)xv7v] = — 15/47r^x^. Rq[x) is used to normalize the interacting correlator 
R(x), i.e. we will examine the ratio /?(x)//?o(.^) for each current. In our lattice simulations, we use 
the free lattice correlators R{){x) to do the normalization to reduce lattice artifacts. 

The point source for quark propagators is put on (0, 0, 0, 0). The boundary condition in the time 
direction is anti-periodic, while in the space direction it is periodic. To avoid different boundary 
condition effects, we fix the time component nt of x in R{x) at zero, i.e. x = {n^^Uy^n^^Q). We take 
all the diagonal points {n,n,n,0) and those lying approximately along the diagonal: {n,n— l,n — 
1,0) and {n,n,n — 1,0). 

Our quenched data set has 40 configurations with lattice size 16^ and gauge coupling j3 =6.1. 
The bare quark mass am^ = 0.015, 0.025 and 0.05. The lattice spacing is 0.08 fm determined from 
the Sommer parameter. The dynamical data set has about 30 configurations for each bare quark 
mass aniq = 0.015, 0.03 or 0.05 and for each topological charge sector \Q\ = or 1. The lattice 
size is 10^ and j8 = 7.2. The lattice spacing is rather coarse: a ~ 0.16 fm. 

3. Results and discussion 

Fig. |l| shows the normalized correlation functions of the three currents from the quenched 
data set for quark mass aniq = 0.015 and 0.05. As we can see, there is a strong attraction in the 
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amq=0.015, quenched -j- 

□ : ,1^, 0* 
o: j", 0" 
X: J^, 1 + 



[] 



: ^ x'S >^ 2< *< I >V *V ^ * 



4 
x/ a 



amq=0.050, quenched 

□ : J^, + 

o: j", 0" 

X: j3, 1+ [|] 

X XX XX *X Y 



4 
x/a 



Figure 1: Comparison of the normalized correlation functions of the three currents from the quenched data 
set. The squares are for which contains a scalar ("good") diquark structure. The diamonds are for 
which contains a pseudoscalar diquark structure. The crosses are for which contains an axial vector 
("bad") diquark structure. 



scalar diquark channel, which is similar to the attraction seen in the pseudoscalar meson channel in 
Refs. [19, 18]. Also we can see a quark mass dependence in this channel. As the quark mass 
decreases, the attraction increases. Since our lattice spacing a is 0.08 fm, the region we are looking 
at is fm < ;c < 0.64 fm. The correlator for has a tendency to go negative at large distance as the 
quark mass decreases. This behavior is very similar to what observed for the correlator of the scalar 
meson in Ref. [|l^]. It was argued in Ref. [|l^] that zero modes are the source of the negativity. We 
will investigate the zero mode contributions in all three channels later. The correlator for is flat 
and stays close to one. This means that there is little correlation among quarks in this channel. 

To see zero mode effects, we separate the 40 configurations in the quenched data set into two 
groups: 35 of them with topological charge Q^O and the other 5 with 2 = 0. i.e. the first group 
contains zero modes while the second one does not. From each group we compute the correlation 
functions. The results are compared in Fig. ^for quark mass aniq = 0.015. The attraction in the 
current at large x is mainly from zero mode contributions. The repulsion in the current j' also 
has a big contribution from zero modes. At quark mass aniq = 0.05, the correlators from the two 
groups agree with each other within error bars (no graph is shown here). It is not surprising that 
zero mode effects become large at small quark masses since the zero mode contribution to the 
quark propagator is proportional to the inverse of the quark mass. For the current 7^, we do not see 
apparent zero mode effects. In Fig. ^ we compare the three correlators obtained from configurations 
without zero modes. They are not so different as in Fig. [T] when there is no zero mode contribution. 
Nevertheless, it seems that there is more attraction in channel and j' than in channel P around 
X ^ Aa = 0.32 fm. 

The results from the two flavor dynamical simulation cannot be compared to the quenched 
results directly because the dynamical simulation is done with fixed topological charge \Q\ = 0, or 
1 . Thus we will only compare the dynamical results from different Q sectors to investigate the zero 
mode effects. 

The comparisons of the correlators from the two topological sectors are given in Fig. ^. There 
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J^, amq=0.015 
□ : Qs^O, K: Q = 



- X 
0? 



(b) 

4 
X /a 




Figure 2: Comparison of correlators from the topological charge 2 7^ sector (squares) and the Q — Q sector 
(fancy squares) in the quenched data set. The graph on the left is for 7^, the one on the right is for j' . 



am„ = 0.015, = 
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Figure 3: Comparison of the correlators of the three currents from the topological charge Q — Q sector in 
the quenched data set. The three channels are not so different when there is no zero mode contribution. It 
seems that there is more attraction in channel and j' than in channel around x ^ Aa = Q32 fm. 



is no difference within error bars between the correlators from the two sectors for both current J 
and . Similarly, no difference is seen for the current J^. In Fig. ^ the 2 7^ sector of the quenched 
simulation contains not only configurations with |2| = 1 but also those with \Q\ > 1. To compare 
with the dynamical simulation results here, we pick out the |2| = 1 and 2 = configurations and 
calculate the correlators again. The results are given in Fig. ^. The difference between the two 
sectors is small, just like the dynamical results in Fig. ^ Comparing Fig. |5| with Fig. ^, we see that 
the attraction in the channel and the repulsion in the channel in the quenched simulation are 
mainly from configurations with big topological charge: \Q\ > 1. 

4. Conclusions 

The quenched simulation shows that the attraction in the channel is the strongest. con- 
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Figure 4: Comparison of correlators from the topological charge 2 = sector (squares) and the |2| = 1 

sector (diamonds) from two flavor dynamical simulations. The graph on the left is for J^, the one on the 
right is for j' . There is no difference within error bars between the correlators from the two sectors. 



J^, amq = 0.015 
□ : |Q|=1, K: Q = 
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Figure 5: Comparison of correlators from the topological charge 2 = sector (fancy squares) and the 
\Q\ = 1 sector (squares) from the quenched data set. The graph on the left is for J^, the one on the right is 
for j' . The difference between the correlators from the two topological sectors is small. 



tains a scalar diquark structure. By comparing the correlators from different topological sectors, we 
found that both the attraction in the scalar diquark channel and the repulsion in the pseudoscalar di- 
quark channel have big contributions from configurations with more than one zero mode. The cor- 
relators for the currents and obtained from configurations without zero modes in the quenched 
data set are similar to each other. Both of them seem to have more attraction around x ~ 4a = 0.32 
fm than the correlator for the current which contains an axial vector diquark structure. The cor- 
relators from the two flavor dynamical data set show no difference between the g = and |Q| = 1 
sector. This confirms that zero mode effects are mainly from configurations with big topological 
charges. In this respect, quenched and full QCD in sectors of low Q are not too different. 

In full QCD, small quark masses suppress high Q configurations. Since they are absent there, 
we suspect that different diquarks in light baryons are not that different. At a minimum, diquark 
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contributions in quenched QCD {not filtered by Q) and full QCD are different, and results from 
quenched simulations may be misleading. 
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